Electron Correlations in a Quantum Dot with Bychkov-Rashba Coupling 



O 

o 

-i— > 
O 

O 



Tapash Chakraborty** and Pekka Pictilalnen*''' 
* Department of Physics and Astronomy, University of Manitoba, Winnipeg, Canada R3T 2N2 
' Departement of Physical Sciences/ Theoretical Physics, P.O. Box 3000, FIN-90014 University of Oulu, Finland 

(February 2, 2008) 

We report on a theoretical approach developed to investigate the influence of Bychkov-Rashba 
interaction on a few interacting electrons confined in a quantum dot. We note that the spin-orbit 
coupling profoundly influences the energy spectrum of interacting electrons in a quantum dot. Inter- 
electron interaction causes level crossings in the ground state and a jump in magnetization. As the 
coupling strength is increased, that jump is shifted to lower magnetic fields. Low-field magnetization 
will therefore provide a direct probe of the spin-orbit coupling strength in a quantum dot. 
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It has long been recognized that a two-dimensional 
electron gas (2DEG) in narrow-gap semiconductors, par- 
ticularly in InAs-based systems with its high values of 
the ^-factor, exhibit zero-field splitting due to the spin- 
orbit (SO) coupling [1]. This coupling is also the driv- 
ing mechanism for making futuristic devices based on 
controlled spin transport, such as a spin transistor [2,3], 
where the electron spins would precess (due to the SO 
coupling) while being transported through the 2DEG 
channel. Tuning of this precession in the proposed spin 
transistor would provide an additional control that is not 
available in conventional devices, but may be crucial for 
the rapidly emerging field of semiconductor spintronics 
[3]. Hence the upsurge of interest in recent years for a 
better understanding of the SO coupling in nanostruc- 
tured systems. 

The spin-orbit interaction in semiconductor het- 
erostructures can be caused by an electric field perpen- 
dicular to the 2DEG. Riding on an electron, this electric 
field will be felt as an effective magnetic field lying in the 
plane of the 2DEG, perpendicular to the 2D wave vector 
k of the electron. We consider the Bychkov-Rashba (BR) 
Hamiltonian [4,5] 
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where E is the confining electric field at the 2DEG, a = 
(<r x , a y ,a z ) denotes the Pauli spin matrices, and c is the 
speed of light. The single-electron Hamiltonian for the 
2DEG with the electric field normal to the interface, E = 
(0, 0, E z ), takes the form 
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where a is the SO coupling parameter which is sample 
dependent and is proportional to the electric field (inter- 
facial and externally applied). Experimentally observed 
values of a lie in the range of 5 - 45 meV nm [1]. The 



energy dispersion then consists of two branches 
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S ± (k) = k 2 ±ak 
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with an energy separation Ago = S + — £ ~ = 2ak for a 
given k. The corresponding wave functions are 
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The spin parts of the wave functions x (k x , k y ) are mu- 
tually orthogonal and (x ± |cr z |x ± ) = 0. Therefore in the 
states X P ± the spins of the electrons lie in the xy-plane 
and point in opposite directions. In addition, 
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i.e., the spins are perpendicular to the momentum 
(k x ,k y ). Spatial alignment of spins therefore depends 
on the wave vector [1]. 

Spin-orbit interaction and electron-electron interac- 
tions are responsible for a variety of interesting effects 
in quantum dots [6]. In this paper, we present a numeri- 
cally exact treatment of the BR Hamiltonian in a system 
of interacting electrons confined in a parabolic quantum 
dot (QD) [7,8] under the influence of an external mag- 
netic field. More specifically, we explore the energy spec- 
tra and magnetization [9] of a few interacting electrons 
in a quantum dot in the presence of SO coupling. It 
should be pointed out that while a large number of theo- 
retical work has been reported as yet in the literature for 
a 2DEG [10] and a QD [11,12] with spin-orbit coupling, 
in most cases, the electron-electron interaction has been 
ignored (or treated within an approximation [13]) due to 
its inherent complexity. 

Let us begin with the single-electron states. Unlike the 
case of a circular quantum dot with hard walls where ex- 
act analytical results for the single-electron energy spec- 
trum (in the presence of the SO interaction) are available, 
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for the more realistic case of a parabolic QD, the en- 
ergy spectrum can only be obtained numerically. In the 
presence of the BR interaction, the Schrodinger equation 
consiste of two parts 
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where -0 is a two-component spinor 
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^ ± iJj-, and v c (r) is the confinement potential. 



We seek a solution of the form 



^ = fr(r)e im i e 



which with Eq. (1) yields, rri| = m^ — 1(= m). In the 
case of a parabolic confinement potential v c = im*<x>Qr 2 , 
the radial equations are 
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where x = r 2 /a 2 , a 2 = h/(m*uj n ), v = e/(2hojo) and 
(3 = m*aa/Ti 2 . When (3 = (i.e., a = 0), Eq. (3) reduces 
to two uncoupled Laguerre equations with solutions 

/ T = f nm = e-^M/^M 



with the energies, 
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(p — e/cA)] z of the SO interaction. The last term gives 
the Zeeman energies of the components of the spinors. 
When (3 = the functions f nm will still be eigenstates 
of the Hamiltonian provided that we replace the single 
particle energies v nm with the expressions 
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where signs a = ±1 correspond to the upper and lower 
components of the spinor. Furthermore, for the angular 
velocity wo related to the harmonic confinement poten- 
tial we have to substitute the effective angular velocity 

1 in 

CI = ojo (l + ojI / (Abjl)) where lu c = eB/(m*c) is the 
cyclotron frequency. 




In the presence of an external magnetic field B the term 
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has to be added to the spinor Hamiltonian Ti. Here the 
first two terms (diagonal) are due to the interaction of the 
orbital motion and the magnetic field. The third (non- 
diagonal) term originates from the vector potential part 
A = -j (— y, x,0) in minimal coupling scheme a/h[& x 



FIG. 1. Energy spectrum of a two-electron InAs quantum 
dot versus the applied magnetic field for different values of 
the Bychkov-Rashba interaction parameter a (meV nm). For 
clarity, at each value of the magnetic field and for a given 
total angular momentum, only the lowest energy is plotted. 
The Zeeman energy is also included. 



For the (3 ^ case which is our main concern here, we 
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use the following expansion 
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For the angular momenta m > we find that the expan- 
sion coefficients satisfy the equations [11] 

iy - «4m)4 =2 *7 + ( n + m + + *? _nc i-i 



and the equations 
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for states with m < 0, and T7+ = l±ea 2 i?/ (Tic). Solutions 
of these eigensystems provide the single-electron energies 
e and the spinor wavefunctions tp [Eq. (2)], which have 
been investigated earlier by several authors in a variety 
of ways [12]. 




FIG. 2. Same as in Fig. 1, but for the InSb quantum dot 
system. 



For a system of interacting electrons we diagonalize 
the many-body Hamiltonian in a basis consisting of non- 
interacting many-body states, which in turn are con- 
structed as antisymmetrized direct products of the two- 
component spinors tp. Since the Coulomb force is inde- 
pendent of the spin orienation we evaluate the sum 

WxMVcouMxM = E (CCI^oullCO. 



A = (77,777), of four terms. Explicit expression for these 
terms for a parabolic QD can be found in Rcf. [8] . 

In our numerical investigations, we choose InAs 
and InSb quantum dots with parameters, m* /m = 
0.042,e = 14.6,5 = -14 and m*/m = 0.014, e = 
17.88,5 = —40 respectively. InSb quantum dots are also 
considered here because of its high g values. In both sys- 
tems, we choose Tiwo = 7.5 meV. The energy spectrum of 
the two-electron state in InAs QD is shown in Fig. 1 for 
various values of the BR coupling parameter a. Similar 
results for InSb QD are presented in Fig. 2. 

The essential feature of the energy spectra at a = is 
that with the increase in magnetic field, the ground state 
moves from J = to J = 2 ( J = m + s z is the total 
angular momentum). This is already well established in 
the literature [9] . This level crossing persists for nonzero 
value of a, but the crossing point shifts to lower magnetic 
fields. This shift of the crossing point can perhaps be 
observed experimentally by a variety of ways, such as 
capacitance spectroscopy, or by transport spectroscopy 
[8]. 

The results for magnetization at the ground state, de- 
fined as M = —dE/dB, where E is the total energy of 
the system, of quantum dots with or without the BR 
interaction is presented in Fig. 3. Magnetization is a 
fundamental thermodynamic quantity that reflects the 
change of the ground state electron energy in a magnetic 
field [14], thereby providing valuable information about 
many-electron dynamics of the QD in a magnetic field. 
We have established earlier that oscillations in magneti- 
zation in a few electron-quantum dot are a direct conse- 
quence of the effects related to the electron-electron in- 
teraction between the two-dimensional electrons confined 
in the dot [9]. A jump in M occurs at a magnetic field 
where the ground state changes from one angular mo- 
mentum to another (Figs. 1, 2). Similar behavior is also 
expected in a nanoscopic quantum ring [15]. With in- 
creasing strength of a, this jump in magnetization at the 
energy-level crossing is pushed to lower magnetic fields. 
For the InAs QD this shift can be as large as ~ 1.5 Tesla 
when a is increased from zero to 40 meV nm. Therefore, 
low-field magnetization measurements of quantum dots 
could be a direct probe of the SO coupling strength. 

In closing, we have developed a theoretical approach 
where the SO interaction is treated via exact diagonaliza- 
tion of the Hamiltonian for interacting electrons confined 
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in a parabolic QD. Coulomb interaction causes energy 
levels to cross and at the crossing point magnetization 
shows a jump. In a magnetic field the strength of the SO 
coupling is proportional to the field (in addition to the 
coupling parameter and the angular momentum). Hence, 
the effect of the coupling is more prominent for slopes of 
the higher angular momenta energy curves. As a conse- 
quence, an increase in the SO coupling strength causes 
the energy level crossings to move to weaker fields and 
the jump in magnetization shows a large shift to weaker 
magnetic fields. This result can be exploited to tune 
the SO coupling strength that might be useful for spin 
transport. Our theoretical approach can be extended to 
include larger number of electrons in the dot. Details will 
be published elsewhere. 
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FIG. 3. Magnetization in the ground state for various val- 
ues of the SO coupling strength and for the two electron QD 
systems. 
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